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We prove sharp I? boundary decay estimates for the eigenfunctions of 
certain second order elliptic operators acting in a bounded region, and of 
their first space derivatives, using only the Hardy inequality. These imply 
L 2 boundary decay properties of the heat kernel and spectral density. We 
deduce bounds on the rate of convergence of the eigenvalues when the region 
is slightly reduced in size. It is remarkable that several of the bounds do not 
involve the space dimension. 



AMS subject classifications: 35P99, 35P20, 47A75, 47B25 
| keywords: boundary decay, Laplacian, Hardy inequality, eigenfunctions, 

Q\ ■ heat kernel, spectral density, spectral convergence. 

1 Introduction 

Let if be a non- negative second order elliptic operator acting in L 2 (U, d N x) subject 
to Dirichlet boundary conditions, where U is a bounded region in H N or even in 
a Riemannian manifold. Let d be a continuous function on U satisfying |Vd| < 1, 
for example the distance from the boundary of U, which may be very irregular. 
We say that H satisfies Hardy's inequality with respect to d if 



jJJl < C * (W) + 



for all / G C^°(U), where Q is the quadratic form of H. The precise value of the 
constant c will be of great importance below, but the size of a is not crucial. 

We are concerned with boundary decay of the eigenfunctions of H, and more 
generally of any functions in the domain of H. Such bounds were first obtained in 
|| |6| by an argument related to that which we use below. The bounds were used 
in to obtain explicit upper bounds on the rate at which 

\X n (U) - X n (U £ )\ 
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vanishes as e — > 0, where X n (S) denotes the n-th Dirichlet eigenvalue of any region 
S and 

U e := {x G U : d(x) > e}. 
In two recent papers Pang |T2L 0| used a different method to obtain a sharp rate 



of convergence of the first eigenvalue as e — > for a Dirichlet Laplacian in any 
simply connected subregion of R 2 , and for certain bounded regions in K N . In this 
paper we return to the method of || and show how to obtain sharp estimates of 
the rate of convergence for all eigenvalues; see Section 4. 

The key is to obtain better boundary decay estimates of the eigenf unctions, in 
terms of 

i/i 2 



L 



' {x:d(x)<e} 

for all / € Dom(if ) and all e > 0, instead of estimating 

l/l 2 



u <P 

for all possible 7 > 0. It is well known to harmonic analysts that the former type of 
estimate is generally sharper than the latter, and we find that it yields the optimal 
power in the subsequent proof of the convergence of the eigenvalues. 

The main theorems of the paper in Section 3 apply to weighted Laplace-Beltrami 
operators acting in regions with irregular and possibly fractal boundaries, but in 
Section 5 we show that the methods can be applied to second order uniformly 
elliptic operators with measurable highest order coefficients. The estimates are 
proved for functions in the domains of the operators, and apply in particular to 
eigenf unctions. In most theorems we prove that we have the optimal power of e in 
the estimates. 

The methods which we use do not require U to be a region in a Riemannian 
manifold. If U is a piecewise manifold obtained by glueing together manifolds of 
the same dimension along certain common edges, the same ideas can be applied 
provided the operator H is defined by means of the appropriate quadratic form. 

In Section 6 we use the results to obtain some new L 2 boundary decay estimates 
for the heat kernel of the operator, and remark that the same methods can be used 
for the spectral density. 

The sharp constant in Hardy's inequality is the only important input to the ar- 
gument, and we refer to || for a recent review of this topic. Here we mention 
only a few outstanding results for H := —A DIR acting in a bounded region U in 
Euclidean space. If U is a simply connected proper subregion of R 2 then Hardy's 
inequality holds with c = 4 and a = by [|, @, Th. 1.5.10]. If U is a convex region 
111 R N then it is a matter of folklore that Hardy's inequality holds with c = 2 and 
a = 0; some proofs are described in ||. Finally, if U has smooth boundary then 
Hardy's inequality holds with c = 2 for some a < 00, H. 
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2 Definitions 



We follow the notation of |J but with somewhat more restrictive conditions on the 
various coefficient functions. Let a be a measurable function on the incomplete 
Riemannian manifold U which is positive almost everywhere and locally L 2 wth 
respect to the Riemannian volume element. Define the weighted space L 2 (U) to 
be the set of (equivalence classes modulo null sets of) functions such that 



/ |/|V<oo. 
Ju 



This and subsequent integrals are evaluated using the Riemannian volume element. 
The introduction of the weight o allows extra applications of our theorems at no 
cost. Let V be a non- negative locally L 1 function on U and let H be the operator 
on L 2 (U) defined formally by 

Hf := -a- 2 V ■ (a 2 Vf) + Vf 

subject to Dirichlet boundary conditions. Rigorously we start from the non- 
negative quadratic form 



W) := / (|V/| 2 + y|/| 2 )a 5 
j u 



which is well-defined on the domain C^°(U) by the hypothesis on a. We assume 
that Q is closable and define H to be the self-adjoint operator on L 2 (U) associated 
with the closure of the form as described in j|, Ch. 4] and pi Section 1.2]. For 



discussions of conditions on a which imply that Q is closable see |L4| and || Section 
1.2]. 

If U is a region in R , a = 1, V = and we choose the Euclidean metric then 
H = —A subject to Dirichlet boundary conditions. Our results are new in this case 
when U is bounded and its boundary dU is irregular, possibly fractal, improving 
on the recent theorems in B |12|, |13|]. 



Our main assumption is formulated in terms of a positive continuous function d 
on U such that \Vd\ < 1, in the weak sense. More precisely, we assume that 

\d(x) - d(y)\ <\x-y\ 

for all x, y G U. This is equivalent to the statement that d has distributional 
derivative Vd G L°° which satisfies |V<i(x)| < 1 almost everywhere in U. One 
might take d(x) to be the distance of x G U from any closed subset S of dU 
or from a closed subset of M\U if U is embedded in some larger Riemannian 
manifold M. We assume throughout the paper that for some constant c > 2 and 
some non-negative constant a the Hardy inequality (HI) 



/mw (g(/)+ 
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is valid for all / G C%°(U), and hence for all / in the domain of the closure of Q. 
Our goal is to obtain a similar but stronger bound for all / G Dom(if ) and hence 
for all eigenfunctions of H. Note that we do not assume that U is bounded or H 
has discrete spectrum. 

A serious difficulty is the fact that we cannot identify the domain of H with any of 
the Sobolev or other spaces. If the coefficients of H or the boundary are irregular 
the operator domain of H changes if we vary o or V within the permitted classes, 
even though the quadratic form domain may be unchanged. The bounds which we 
obtain in Theorems 4 and 7 bear some relationship with Morrey space estimates, 
already known to be of great importance in the theory of elliptic operators, [0, [H], 



3 The main theorems 

Our estimates involve a positive parameter s, and various other constants which 
depend only on c > 2, in a way which we make explicit. Given e > we put 

u(x) := (max{d(x), e})^ 1 ^ 

for all x G U . 

Lemma 1 If f G Dom(i7) and s > then 

{Hf^ 2 f) + s\\uf\\l\<c^\\{H + s)f\\ 2 \\{H + a) 1 l c f\\ 2 . 

Proof Using HI and j|, Lemma 4.20] we have 

^ 4 <(cT 2 ) 2/c <{c 2 (# + a)} 2/c 



so 

and 

Hence 



< (H + a)-^ c u\H + a)-^ c < c^ c I 
= ((H + s)f,u; 2 f) 

= ((H + s)f,u 2 (H + a)- 1 / c .(H + a) 1 / c f) 
<<?' e \\(H + 8 )f\\ 2 \\(H + a) 1 '°f\\ a . 
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Lemma 2 If f G Dom(Q) and // G W 1 ' ^?/) i/ien /// G Dom(Q) and 

0(a*/)<2||a*||2oW) + 2||Vai||2oI|/||^ 

Proof If / G C™{U) then fjtf G W,, 1 ' 00 ^) C Dom(Q) and 

QM) = ^(|/iV/ + /V//| 2 + Vl/| 2 )^ 2 

(2 / , 2 |V/| 2 + 2|/| 2 |V/i| 2 + V|/| 2 ) * 2 



< 

< 2|| / u|| 2 <) Q(/) + 2||v^|| 2 0O |mi 2 . 

If / G Dom(Q) then the fact that fif G Dom(Q) and the validity of the same 
estimate both follow from the lower semi-continuity of Q. 

Lemma 3 If f G Dom(if ) then 

lu ~ C2/1{H + a)fU{H + a)1/Cfh + i |V " |2|/I V " 

Proof By HI and Lemma 2 we see that u;/ G Dom(Q) and 

^ 2 |/|V 



c 2 rf 2 



<Q(u;/) + a||u;/|| 2 . 



Secondly 



Q(ujf)-~(Hf,u; 2 f)-~(u; 2 f,Hf) 



jT {|v(^/)| 2 - ±v/ ■ v(w>/) - iv(^ 2 /) ■ V/} a 2 

/ |Vc| 2 |/|V. 



Hence 



Q(uf) + a\\uf\\ 2 
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= 1 (<#/, u 2 f) + a||^/|| 2 ) + i ((u 2 f, Hf) + a|K|| 2 ) + jf | V^| 2 |/| V. 

The proof is completed by combining the above two formulae with the bound of 
Lemma 1. 

For some comments on the optimality of the estimates in the following theorem 
see Example |5|, Example [| and the note after Corollary ||. 

Theorem 4 If f G Dom(if ) then assuming HI we have 



{x:d{x)<e} d? 



ltl \ 2 <coe 2 / c \\(H + a)f\\ 2 \\(II + a)V c f\\ 2 



for all e > 0, where Cq := c 2+2 / c . Hence 

[ \f\ 2 a 2 <c,e 2+2 l c \\{H + a)f\\ 2 \\{H + a)^f\\ 

J{x:d(x)<e} 

for all e > 0. 

Proof We rewrite Lemma 3 in the form 

y|/|V 2 <c 2 / c ||(i/ + a)/|| 2 ||(// + a) 1 /7l|2 



u 

where 

Y ■ = 

c 2 d? 



If d(x) > e then 



iVcul 2 < c- 2 d- 2 - 2 ' c ^ 



c 2 d 2 

so Y{x) > 0. On the other hand if d(x) < e then 

y = ^> 1 



c 2 d 2 ~ c 2 e 2 l c d 2 ' 
Thus 

/ \IL a 2 < C V/ C / F|/|V 

^{x:d(a;)<e} GT J[7 

< c 2+2 ^ c \\(H + a)f\\ 2 \\(H + a)^f\\ 2 . 
The second statement of the theorem is an immediate consequence of the first. 
Example 5 If c = 2 and a = 1 then the theorem states that 

/ \f\ 2 <8s 3 \\(H + a)f\\ 2 \\(H + a)^fh 

J{x:d(x)<e} 

for all / G Dom(if), which is the L 2 analogue of f(x) = 0(d(x)) as d(x) — > 0. In 
particular suppose that U C R w is bounded with a smooth boundary <9?7, and let 
/ be a generic function in C^°(U) which vanishes on dU. Then / G Dom(if) and 
the power 3 of e above is optimal. 

Example 6 Let U : = (0, oo), d(x) := x and cr(x) := x a ^ 2 where < a < 1. Then 
the operator if is given formally by 

Hf(x) := 



subject to Dirichlet boundary conditions at x = 0. The quadratic form Q has 
domain W 1,2 ((0, oo), x a dx). A standard result, [0, p. 104], states that the strong 
Hardy inequality holds with c = 2/(1 — a) 
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Now let / be a smooth function on (0, oo) which vanishes for x > 2 and equals 
x i-a f or < a < 1. It is easy to prove that / G Dom(if) C W 1,2 ((0, oo),x a dx). If 
< £ < 1 then one also has 

J{x:d(x)<e} d 2 3 — a 

Therefore the power 2 + 2/c in Theorem |] is optimal. 

Corollary 7 If g(s) is a monotonically decreasing C 1 function on (0, 5} which 
vanishes for s = 5 then 

I g(d)\f\ 2 o- 2 < co f \g'{s)\s 2+2 /^s.\\{H + a)f\\ 2 \\{H + a) l / c f\\ 2 (2) 
Ju Jo 

for all f G Dom(H), provided the integral on the RHS is finite. 
Proof We have 

g(d)= [ 5 \g'(s)\ds 
Jd 

for all d G (0, 5}. The corollary follows by applying Fubini's theorem to 

£ {Xx«.| 9 'MII/|V}d s 

where x stands for the characteristic function of a set. 

Note Let H := — £±dir acting in L 2 (U,dx) with a := 1, where U is a bounded 
region in R^, and let d be the distance to the boundary dU. If g(s) = o(s~ 2 ~ 2 ^ c ) 
as s — ► then (|2|) is equivalent to 

/ g(d)\f\ 2 <c (2 + 2/c) f g{s)s l+2 / c ds.\\{H + a)f\\ 2 \\{H + a) 1 / c f\\ 2 (3) 
Ju Jo 

which may be compared with the pointwise bound 

\f(x)\ = 0{d{x)W e ) (4) 

as d{x) — > 0. If dU is C 2 then c = 2 and this pointwise bound is sharp for the 
first eigenfunction of H . However, no such pointwise bound exists for arbitrary 
functions in the domain of H . Moreover, if dU is fractal it is not clear that 
is the correct pointwise analogue of (|3|), nor indeed that there is any pointwise 
analogue. 

Our next task is to obtain comparable estimates for |V/|. This necessitates intro- 
ducing the continuous function on r : U — > [0, oo) defined by 

f e- x l c if < d(x) < e 

t{x) := | ^^-^((l + c)e - d{x)) if e < d{x) < (1 + c)e 
[ otherwise. 
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It is immediate that 

< r < u < e~ 1/c 

|Vr| < c^e- 1 - 1 ' 

supp(r) C {x : < d < (1 + c)e}. 



Theorem 8 If f £ Dom(H) then assuming HI we have 

f \Vf\ 2 a 2 <c l e^\\(H + a)fh\\(H + a)^f\\ 2 . 

J{x:d(x)<e} 

for all e > 0, where 

Cx : = c 2 ' c + c 2 / c (l + c) 2+2 / c . 

Proof We have 

e- 2/c f |V/|V 

J{x:d(x)<e} 

= f |V(r/)|V 
< Q(r/) 

where r/ G Dom(Q) by Lemma 2. By the same argument as in (p]) of Lemma 3 
this equals 

21 f|2J 



1 -(Hf,T 2 f) + ±(T 2 f,Hf) + /jVrH/lV 
< IIF/Hall^CF + o)- 1 / 6 !! ||(iy + a) 1 / c /|| 2 + /' |Vr| 2 |/| 



,2 a 2 

< c 2 /1F/|| 2 ||(F + a)^f\\ 2 + c- 2 e- 2 - 2 ^ I \f\ 2 a 2 

J{x:d(x)<(l+c)e} 

< dI|( J ff + a)/|| 2 ||( J H'H-a) 1 /Vl| a 
using Theorem 4. 

Note By extending the calculation of Example ^| one sees that the power of e in 
the above theorem is optimal. The choice of r in the proof is certainly not optimal, 
so neither is the value of C\ obtained. 

Corollary 9 If Hf = Xf and \\f\\ 2 = 1 then 

[ |/|V<c e 2+2/c (A + a) 1+1 / c (5) 

and 

[ |V/|V< Cl e 2 / c (A + a) 1+1 / c 

Jix:d(x)<e\ 



for all e > 0. 
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Proof These follow directly from Theorems 4 and 6. 

Note If we insert the eigenfunction / directly into HI we obtain 

/ |/|V<cV(A + a) 

J{x:d{x)<e} 

which is exactly what is obtained by interpolating between fl5|) and the trivial 
estimate 

/ |/|V < l. 

J{x:d(x)<e} 

This supports the conjecture that the constant cq in Theorem |4] is optimal. 

Corollary 10 If H := —Adir in L 2 (U,d 2 x) where U is a simply connected proper 
subregion o/R 2 and 

d(x) := dist(x, dU) 

then 

[ \f\ 2 <32e^\\Hf\\ 2 \\H 1 ^f\\ 2 

J{x:d(x)<e} 

and 

[ |Vjf < lUe l ' 2 \\Hf\\ 2 \\H^f\\ 2 

J{x:d(x)<e} 

for all f G Bom(H). 

Proof We may put c = 4, a = and a = 1 in Theorems 4 and 6 by Q, |], Th. 
1.5.10]. 

4 Perturbation of the domain 

In this section we use the results above to consider the effect on the spectrum of 
H of replacing the region U by a slightly smaller region U £ such that 

{x e U : d(x) >e}CU £ CU. 

If \ n (U e ) denote the eigenvalues of the operator H e defined by restricting H to 
L 2 (U e ) where we again impose Dirichlet boundary conditions, then variational 
arguments imply that 

Xn(U) < X n (U £ ) 

for all n and e > 0, and our goal is to find quantitative estimates of the difference. 
The constants c n below all depend only on a, c, Co, c\ and n. 

Let \i : U — > [0, oo) be defined by 

[ if < d(x) < e 

H{x) := < (d(x) —e)/e if e < d(x) < 2e 
[ 1 otherwise. 

so that < \l < 1, |V/i| < e^ 1 and [l has support in U e . 
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Lemma 11 There exists a constant c 2 > such that if f G Dom(H) then 
Qfrf) < Q(f) + e 2/c c 2 \\(H + a)fh\\(H + a) l / c f\\ 2 . 

Proof Putting S := {x : e < d(x) < 2e} we have 

Q(»f)-Q(f) < I |V(/x/)|V 
Js 



< 2 / /i 2 |V/|V 2 + 2 f |V/i| 2 |/|V 

Js Js 

< l\ |V/|V + 25- 2 / |/| 



< ^calKff + ^/llalKff + o) 1 /"/!!,, 

by Theorems 4 and 6. 

It is crucial to the application of our next lemma that < e 1+1 ^ c < e 2 l c provided 
< e < 1, so the error is actually smaller than that of Lemma 9 as e — ► 0. 

Lemma 12 There exists a constant C3 > such that if f 6 Dom(ff) then 

\\fh> \Wh> \\fh-c 3 e 1+1 / c \\(H + a)f\\l /2 \\{H + afl c f\\y 2 

Proof The first inequality is elementary. We also have 



/ (i-/^) 2 l/|V 

Ju 



( , -^ 2 '- 2 ' 2 
< I mv 

J {x:d(x)<2e} 

= e^cl\\(H + a)f\\ 2 \\(H + a)^f\\ 2 
by Theorem 4. The second inequality of the lemma follows. 

The case n = 1 of the following theorem with the sharp power e 1 ^ 2 corresponding 
to c = 4 was already proved for proper simply connected subregions of R 2 in ||13|| , 



by an entirely different method which seems not to extend to higher eigenvalues. 
Theorem 13 There exist constants c n for all positive integers n such that 

Xn(U)<X n (U £ )<X n (U)+c n e 2/c . 

Proof this follows || Th. 22] closely. 

5 Elliptic operators 

In this section we extend the earlier results to second order uniformly elliptic 
operators in divergence form with possibly measurable second order coefficients, 
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making use only of the ellipticity constant of the operator. Throughout the section 
we put o"=l and integrate with respect to Lebesgue measure. 

Let U be a bounded region in with C 2 boundary and let 

d(x) := dist(x, dU) 

so that 

for some a > and all / G Wq ,2 (U). Now let 

df\ 



Hf(x) :=-E A (««(*) 



subject to Dirichlet boundary conditions in L 2 (U), where 

1 < a(x) < a 2 

for all x G U, and we interpret H as a self-adjoint operator using the theory of 
quadratic forms as usual. If we put d(x) := a~ 1 d(x) then 

^ ^ dd dd 



for all rr G £/ and 

/^<4« 2 (W) + a||/|| 

for all / G Wo'^C/), where Q is the quadratic form asociated with H. 

We are now in a position to apply the theory of the paper to the pair H, d with 
c := 2a. 

Theorem 14 There exists a constant cq such that if f G Dom(H) then 

I „ J/| 2 <Co£ 2+1 / a ||(// + a)/|| 2 ||(// + a) 1 /^)/|| 2 . ( 6 ) 

./{x:d(x)<e} 

Theorem 15 There exists a constant c\ such that if f G Dom(if) £/ien 

/ ^ jV/| 2 <c l£ 1 / 1 |(/ / + a)/ || 2 ||( // + a) V(2a) / || 2 . (7) 

We next suppose that U £ is a region satisfying the same conditions as in Theorem 
11, and define X n (U £ ) in a similar manner. 

Theorem 16 There exist constants c n for all positive integers n such that 

X n (U)<X n (U £ )<X n (U)+c n e 1 / a . 
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In each case we conjecture that the power of e is optimal. The three theorems 
can be proved in two ways. We may adapt the proofs of this paper, replacing the 
weighted Laplacian by a more general second order elliptic operator. Alternatively, 
we may apply the theorems of the paper, but using a Riemannian metric and weight 
adapted to the choice of the second order coefficients, as described in . Namely 
if gi j(x) is the matrix inverse to a* J (x) then the Riemannian metric 

9i ,j{x)dx' l dx : > 

is Lipschitz equivalent to the Euclidean metric in U. Indeed the Riemannian 
distance function is bounded between or x d and d. If also 



a(x) := det (gij(x)) ' 

then 

/ l/l'V = / \l\'A"x 
Ju Ju 

where the integrals on the left are with respect to the Riemannian measure and 
d N x is the Lebesgue measure. Hence 

Hf := -a 2 V-(a 2 Vf) 




We finally remark that the above theorems can be localised. Suppose H is sub- 
elliptic but the hypothesis 1 < a(x) < a 2 holds for x such that dist(x, S) < j3, 
where (3 > and S is some closed subset of dU. We only assume that dU is C 2 in 
the /3-neighbourhood of S. If we put 

d(x) = a" 1 min{dist(a;, S),(3} 

then 

^ , dd dd 

for all x e U, because the gradient of d vanishes outside the /3-neighbourhood of 
S. The proof of 

/^<4« 2 (Q(f) + a\\f\\t) 

for all / G W ' 2 (U) involves the same arguments as in 0, concentrating on the 
region {x e U : dist(x, S) < (3}. 
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6 Heat kernel and related bounds 



If K(t,x,y) is the heat kernel of a uniformly elliptic second order operator H 
written in divergence form and acting in L 2 (U, d N x) subject to Dirichlet boundary 
conditions, then it is known that 



< K(t,x,y) < c 2 t 



-N/1 



for all x,y G U and all t > 0. If H = —Adir then we may even take c 2 = (4ir)~ N / 2 ; 
see || for an account of the relevant heat kernel bounds. We are interested in L 2 
boundary decay properties of the heat kernel and spectral density which bear some 
relationship with the 'intrinsically ultracontractive' pointwise bounds obtained un- 
der much stronger assumptions and with much less control on the constants in B 
Chapter 4]. Throughout this section we assume (Q), (^) and (H); the constants 
in our bounds depend on these constants and on N in a manner which is easy to 
make explicit. 

Theorem 17 Under the assumptions ffi) and ^) we have 

K{t, x, y) 2 d N x < c 3 e at {e 2 /t) 1+1 / {2a k- N / 2 (9) 



' {x:d(x)<e} 

for all y G U and all t > 0. If also U is bounded then 

[ K(t, x, x)d N x < c A {e 2 /t) l+l ^ 2a h- N ' 2 (10) 

J{x:d(x)<e} 

for all e > and all < t < 1. 

Proof Denoting the left-hand side of (pf) by / we have 



/=/ \e- m 5 y \ 2 d N x= \e- Ht ' 2 g\ 2 d N 



y | 

{x:d(x)<e} J{x:d(x)<e} 



X 



where 5 y is the delta function at y G U and g := e Ht / 2 5 y . Using the semigroup 
property we have 



\\g\\l = ^ K(t/2, x, y) 2 d N x = K{t, y, y) < c 2 t~ 
Applying (|6|) and then the spectral theorem we deduce 



■N/2 



I < cvE 2+1 l a \\{H + a)e- Ht / 2 g\\ 2 \\{H + ayl^e- Ht l 2 g\\ 2 
: 2+1/a e at \\(H + a)e- {H 

< C3£ 2 + 1 /« e ^ t -l-l/(2a)-iV/2 



< c e 2+1 ^e at \\(H + a)e^ H+a ^ 2 \\ || (H + a ) 1 ^ 2a) e~ <yH+a)t l 2 \\ \\g\\ 2 2 



We adopt an alternative strategy to prove the second inequality. Let {\ n }^ =l be the 
eigenvalues of H written in increasing order and repeated according to multiplicity, 
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and let {0„}^L 1 be the corresponding normalised eigenfunctions. It is known that 
there exist positive constants and a 2 such that 

am 2 / N <\ n < a 2 n 2 ' N 

for all n. Also 



K(t,x,x) = Y. e 



n=l 

for all a; e U and t > 0. Denoting the left-hand side of ([[(]) by J we deduce that 
for < t < 1 

j = e e ~ Xnt / i^rd^x 

n= l J{x:d(a;)<e} 
00 

< £e- A "*c £ 2+1 / a (A n + a) 1+1 /( 2a ) 

n=l 

00 

< c Q e 2+1 ^ a E e-^^^^asn 2 ^ + a) 1+1 /( 2a ) 

< C 4 £ 2+1 / a t- 



n=l 

,2+l/aj.-l-l/(2a)-jV/2 



where for the last line we compared the sum with the corresponding integral. 

The estimate (|1(]) is not asymptotically optimal as e, t — > 0, even for H := — A^/^ 
acting in L 2 (0, 00). In this case we have 

K(t,x,x) = (47it)- 1/2 (l-e- x2/t ) 
- (4vr)- 1 /Vt- 3 / 2 

if < x 2 « t by the reflection principle || pl07]. Therefore 

K(t,x,x)dx ~ (3Qtt)- 1/2 sH- 3/2 

Jo 

if < e 2 << t. However Theorem 15 with a = 1 only yields 

/ K(t, x, x)dx < c 4 e 3 t~ 2 . 
Jo 

The following is a possible reason for this failure. Optimal estimates on eigenfunc- 
tions associated with highly degenerate eigenvalues can be much worse than one 
expects for typical eigenvalues. Our proof uses a bound for every eigenfunction 
which takes no account of this fact, so when summed up it is not surprising that 
the resulting heat kernel bound is not optimal. 

We note that the same method may be used to obtain upper bounds on 

„ 00 

J{x:d{x)<e} n=1 
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Another approach may be used to obtain upper bounds on quantities associated 
with the spectral density. Let E\ be the spectral projection of H associated with 
the interval (—00, A), and let e(X,x,y) be its integral kernel. Then 



N(X) = / e(X,x,x)dx 
J u 

where N(X) is the number of eigenvalues of H less than A. If e > we put 



N(e,X):= / e(A, x, x)dx. 

J{x:d(x)<e} 



Theorem 18 Under the assumption (j^) we have 

N(e, A) < c e 2+1/a (X + a) 1+1/(2a) A(A) 
for all A > and e > 0. 
Proof If / = E\f then (|6]) implies 

/ \f\ 2 <c e 2+1 / a (\ + ar^\\f\\l 

J {x:d{x)<e} 

which is equivalent to the operator inequality 

WQeExW 2 <c e 2+1 / a (X + a) 1+1 ^ 

where Q e is the projection given by multiplying by the characteristic function of 
{x : d(x) < e}. We have 

N(e,X) = ti[Q £ E x Q £ ] 

= tr [(Q £ E\)E\(E\Q £ )] 

< ||Q eJ E7 A || 2 tr[JB? A ] 

< c e 2+1/a (A + a) 1+1 /( 2a )iV(A). 

We finally comment that lower bounds on integrals associated with the spectral 
density have recently been obtained by Safarov, using a coherent state method, 
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